Flow of wet granular materials 
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The transition from frictional to lubricated flow of a dense suspension of non-Brownian particles 
is studied. The pertinent parameter characterizing this transition is the Leighton number Le = ^ 
which represents the ratio of lubrication to frictional forces. The Leighton number Le defines a 
critical shear rate below which no steady fiow without localization exists. In the frictional regime 
the shear flow is localized. The lubricated regime is not simply viscous: the ratio of shear to normal 
stresses remains constant, as in the frictional regime; moreover the velocity profile has a single 
universal form in both frictional and lubricated regimes. Finally, a discrepancy between local and 
global measurements of viscosity is identified, which suggests inhomogeneity of the material under 
flow. 

PACS numbers: 83.80.Hj,83.60.-a,47.55.Kf 



In recent years, the flow behavior of granular matter 
has been the subject of considerable controversy. Simple 
questions such as whether viscosity can be properly de- 
fined forgranular systems are still the subject of hot de- 
bate P,U A notion of viscosity is necessary for many ap- 
plications to predict the resistance to flow. Perhaps even 
more crucial, in view of its importance in geophysics and 
civil engineering, is the resistance to flow of wet granu- 
lar materials. Nevertheless, the number of studies on wet 
granular matter is almost negligible compared to that for 
sand. The usual picture is that for an interstitial fluid at 
low viscosity, the material behaves similarly to dry granu- 
lar systems (with frictional contacts between the grains), 
and for an interstitial fluid at high viscosity the material 
behaves similarly to viscously dominated systems (with 
lubricated contacts between the grains) H, ^ IE IS - 

If this picture were true, it is of importance to investi- 
gate what parameters determine the frictional to viscous 
transition, and whether a viscosity can be defined. 

The purpose of this Letter is to answer some of these 
questions. We find that the transition between frictional 
and lubricated flow regimes is characterized by a viscos- 
ity bifurcation 0, |1|, and is governed by the Leighton 
number Le — with rjs the interstitial fluid viscos- 
ity, 7 the shear rate and a the total shear stress. For a 
given interstitial fluid, Le fixes the critical shear rate jc 
at which the transition between the two regimes occurs. 
In the frictional regime, no steady flow exists without lo- 
calization. We also find that the viscous regime is similar 
to the granular flow regime posing the same problem for 
the definition of the viscosity as for dry granular matter 

We study the rheological behavior of a paste (a dense 
suspension) composed of non-Brownian monodisperse 
spherical particles immersed in a Newtonian fluid of the 
same density at a volume fraction (j) of 60%. The par- 
ticles are monodisperse spherical polystyrene beads (di- 



ameter 0.29 ± 0.03 mm, density 1.04 g.cm~^). To avoid 
sedimentation or creaming effects we match the density 
of the interstitial fluid with the density of the particles. 
Adequately salted water {rjs = 1 mPa.s) and Rhodorsil 
10646 silicone oil (77^ — 2300 mPa.s) are perfectly den- 
sity matched. The oil is miscible with water: we mix 
it with salted water in order to have an isodense liquid, 
whose viscosity varies between 1 and 2300 mPa.s. These 
experiments are completed with pastes prepared with sil- 
icone oils of variable viscosities ^] . The rheology is done 
with a vane-in-cup geometry on a commercial rheometer 
that imposes either stress or shear rate. The vane geom- 
etry we use is equivalent to a cylinder (diameter 16 mm, 
height 52 mm) with a rough lateral surface on the scale 
of the beads, reducing the slipping of granular materials 
IS . For the same reason, the inside of the cup (diame- 
ter 26 mm) is covered with the granular particles using 
double-sided adhesive tape. 

In the first set of experiments, we use a 20 mPa.s oil 
as the interstitial fluid. Fig. ^ shows the evolution of 
viscosity with time, for different applied stresses. We 
find that for stresses smaller than a critical stress Uc = 
5 ± 3 Pa, the viscosity of the sample increases in time 
until the flow is halted altogether, which corresponds to 
an inifinite steady-state viscosity For a stress only 
slightly above ctc, the viscosity decreases with time; for 
long times it reaches a (low) steady state value r]o- This 
implies that at the critical stress, the steady-state viscos- 
ity jumps from an infinite to a finite and low value at cTc- 
This behavior implies the existence of a bifurcation of 
the viscosity, which abruptly passes from a flowing state 
above CTc to a jammed state below CTc, as observed pre- 
viously for yield stress fluids and dry granular materials 
0,I1|- Apart from a critical stress, the bifurcation also 
identifies a critical shear rate jc = 0.4 ± 0.1 s~^, beyond 
which steady flows are possible under an imposed stress. 
If shear rate rather than the stress is imposed, below jc 
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FIG. 1: (a) Viscosity bifurcation : under an imposed stress 
the viscosity either grows in time, or decreases; therefore the 
steady state viscosity jumps to infinity at a critical stress. 
This allows us to define both the critical stress and the critical 
shear rate. Before each experiment, the material is presheared 
during 30 s at 30 to obtain a reproducible initial state, (b) 
Flow curve (shear stress vs. shear rate) for 20 mPa.s silicone 
oil, both at an imposed macroscopic shear stress and shear 
rate. The shadow zone is the statistical error bar. 
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FIG. 2: Critical shear rate vs. interstitial fluid viscosity, for 
polystyrene beads in a Newtonian liquid. The line is a power- 
law fit to the isodensity values with a slope of -0.99 ± 0.10. 



the measured shear stress is almost independent of the 
shear rate, which is the hallmark of quasistatic granular 
(frictional) flow P, |^. For high shear rates (7 > 7c), 
stress increases linearly with increasing shear rate, as for 
a viscous fluid (Fig. ^). These observations allow us to 
identify the transition between two different flow regimes. 

To find out what determines the transition between 
the two regimes, we vary the viscosity of the intersti- 
tial fluid. As before, we determine the critical stress and 
critical shear rate are determined from the viscosity bi- 
furcation. Fig. 13 summarizes the main results. The 
critical shear rate 7c is inversely proportional to the vis- 
cosity of the interstitial fluid, whereas the critical shear 
stress CTc turns out to be constant (within experimental 
uncertainty). The former result is consistent with that 
of Prasad and Kytomaa JJ and Ancey and Coussot Q , 
who only used two fluids of different viscosities. 

We conclude from these data that the pertinent param- 
eter characterizing the transition is the Leighton num- 
ber Le — , which represents the ratio of lubrication 
to frictional forces. Ancey and Coussot previously re- 



ported the inverse proportionality with respect to the 
frictional forces (SjjlQ]- By varying the lubrication forces, 
we conclude that the transition is entirely characterized 
by the Leighton number. The critical stress CTc does not 
vary signiflcantly (ctc = 5 ± 3 Pa), so that Lcc has the 
same order of magnitude for tJs varying from 10^'^ to 2.3 
Pa.s {Lec w (7 ± 5) 10~^). The order of magnitude of 
the critical stress is roughly that of the low shear vis- 
cosity of the suspension multiplied by the critical shear 
rate. The Krieger-Dougherty model [ll| for hard spheres: 

77 = (1 - (j)/(j)m)~^'^'''"", gives crc= riic « 1 Pa. This 
gives the correct order of magnitude of (Tc but also ex- 
plains why the critical stress remains constant: this fol- 
lows from combining 7c oc I/77S and ?/ oc 77^ [T^ . 

All of the rheometric experiments agree with the hy- 
pothesis that the low-shear, low interstitial fluid viscosity 
regime is frictional, in the sense that the material be- 
haves similarly to dry granular matter. The remaining 
question is whether the second regime is viscous. Surpris- 
ingly, measurements of the first normal stress difference 
A'^i (neglecting the second normal stress difference) using 
a plate-plate geometry (diameter 40 mm) show that the 
normal and the viscous stresses are proportional in both 
flow regimes. For the plate-plate geometry, the critical 
shear rate is found to be very similar to what found in 
the Couette ceU (7c = 0.5 ± 0.2 s^^). Fig. shows 
that the ratio a/Ni does not vary significantly with the 
shear rate a/Ni — 0.39 ± 0.15 in both regimes, and over 
4 decades in shear rate. Previous results in the 'viscous' 
regime are consistent with our findings 0, 0, ^3 j despite 
the fact that others B have reported different results. 
Brady and Morris have shown theoretically that a 
shear stress and a pressure proportional to shear rate 
may result from hydrodynamic interactions between the 
particles; however in order to generate normal stress dif- 
ferences, hard sphere contacts must be incorporated. As 
a consequence, there is no reason for the ratio a/Ni to be 
equal in the frictional and lubricated regime as observed 
experimentally. This surprising observation is reinforced 
by the finding that the ratio a/Ni is equal to the inter- 
nal static friction coefficient of the dry granular material 
= 0.40 ± 0.04 (as obtained from the slope angle of a 
heap of dry beads with /is = tan 6*). We conclude that, 
in the lubricated regime, the material has a number of 
characteristics of dry granular materials, while it simul- 
taneously dissipates viscously in the interstitial fiuid (as 
follows from 77 cx rjs)- 

These puzzling observations raise the following ques- 
tions: what happens in the flow? What is the distinction 
between the two regimes? To address these questions, we 
carried out experiments with a velocity controlled "MRI- 
rheometer" , which allows for a direct measurement of the 
local velocity distribution in a Couette geometry [l5l |. 

The main results from these MRI measurements are 
that the velocity proflles are roughly exponential as in 
dry granular materials , and that they occupy only a 
small fraction of the gap at low rotation rates, i.e. in the 
frictional regime: we observe shear localization. How- 
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FIG. 3: Shear/normal stresses ratio a/Ni as a function of the 
shear rate for polystyrene beads in 20 mPa.s sihcone oil (sohd 
volume fraction: 58%), at imposed shear rate. Both plates 
are covered with sand paper to avoid wall slip. The shadow 
zone is the statistical error bar. The horizontal line appears 
off center because of the logarithmic scale. 
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FIG. 4: Velocity rescaled with the velocity of the inner cylin- 
der (inner cylinder radius Ri — 4.15 cm, outer cylinder radius 
Re = 6 cm, height 11 cm; we again use the 20 mPa.s Rhodor- 
sil silicone oil). For technical reasons, Vi can be varied either 
between 0.01 and 3.91 cm/s, or between 0.43 and 43.5 cm/s. 
As in the macroscopic rheology experiments, we preshear the 
material at a rotational speed Vi=43.5 cm/s during 30 s. For 
the low velocity setup, we preshear the material at the maxi- 
mum possible rotational speed Vi, i.e. 3.91 cm/s. The rough 
inner cylinder is driven at a rotational velocity Vi ranging 
between 0.013 and 43.5 cm/s, yielding overall shear rates be- 
tween 0.006 s~^ and 20.9 s~^, the critical shear rate 7c being 
0.4 ± 0.1 s"^ (or Vc = 1.04 ± 0.26 cm/s). 



ever, if the rotation frequency is increased, a surprising 
behavior is observed: contrary to what happens for dry 
granular matter, the higher the rotation rate, the larger 
the fraction of the paste that is sheared (Fig. 0]). In the 
lubricated regime, beyond jc = 0.4 ± 0.1 s~^, the whole 
sample is sheared (7c measured with the MRI is the same 
as found in the rheology) [T^ . 

In the frictional regime, the behavior is only different 
in the sense that a smaller fraction of the material is 
sheared. In this regime the reduced velocity V{R)/Vi 
[Vi being the velocity of the rotating inner cylinder) can 
be collapsed onto the same universal curve when plotted 
as a function of the rescaled coordinate {R — Ri)/dc{Vi) 
(Fig. 13). The length dc{Vi) simply gives the extent of 
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FIG. 5: Velocity profiles measured in a Couette geometry. 
We plot V/Vi vs. {R — Ri)/dc{Vi) where Ri is the inner cylin- 
der radius and dc an adjustable parameter; we also plot two 
velocity profiles extracted from ^18| for the shearing of a dry 
granular material (1 mm mustard grains) in a Couette geom- 
etry (Ri — 3 cm. Re = 6 cm). The line is an exponential fit: 



V/Vi = e ' "^cf^i' ; the dotted line is for a power-law fluid fit 
with exponent n = 0.13. Inset: dc vs. Vi for a 3.91 cm/s 
preshear (open triangles) and a 43.5 cm/s preshear (squares). 



the material that is sheared. In addition we find that 
this universal rescaling of the velocity profile also applies 
to dry granular materials: two velocity profiles from ^ 
and [18] fall along the same universal curve (Fig. [S)), 
underlining the universality of the roughly exponential 
velocity profiles which goes beyond granular pastes. The 
inset to Fig. [S] gives the extent of the sheared region dc 
as a function of the rotation rate. Starting from low Vi, 
dr. increases and eventually fills the whole gap for Vi > Vc 



This last observation provides a natural explanation 
for the macroscopic rheology data. In the first, fric- 
tional regime, no steady flow without localization can be 
achieved at imposed shear rates : there is a coexistence 
between a sheared and an unsheared region, which results 
in a constant stress, not unlike stress plateaux and the 
corresponding shear bands observed for certain surfac- 
tant and polymer solutions . When the sheared region 
has invaded the gap, there is no longer coexistence and 
the stress increases again with increasing shear rate. The 
transition between the frictional and lubricated regimes 
then happens at the end of the coexistence. 

The interesting question is then whether a viscos- 
ity can be defined for the material, in the lubricated 
regime at least, that completely characterizes its resis- 
tance to fiow. If we suppose our material is a contin- 
uum medium in stationary fiow, momentum conservation 
leads to cr{R) = (TiRf/R^ (where ai is the total shear 
stress on the inner cylinder), independently of the con- 
stitutive equation of the material: the stress varies within 
the gap. The shear rate is the spatial derivative of the 
velocity profile. Therefore it varies within the gap and a 
single MRI experiment allows to recalculate a complete 
stress-shear rate curve. The result is shown in Fig. |^ 
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FIG. 6: Stress vs. shear rate calculated from the velocity 
profiles for different rotational velocities and measured macro- 
scopically (same data as Fig. ^). 



for different rotation velocities. There are two surprises: 
first, the data are not consistent between different MRf 
experiments, and second, they are inconsistent with the 
macroscopic rheology data |j21j. We are therefore forced 
to conclude that the velocity profiles are not consistent 
with the macroscopic rheometric data, both when there 
is an unsheared region in the material and when the ma- 
terial is entirely sheared. The possible reasons for this 
is that either the material is inhomogeneous, or that no 
simple constitutive equation relating shear stress to shear 
rate exists for the material (as in dry granular materials 

The question is then: what happens in the lubricated 
regime? We can reasonably assume that, even if force 
chains are present in the sheared system, momentum con- 
servation is likely to hold, at least on average, therefore 
a oc (Ti/R^. If we combine this with the roughly expo- 
nential decay of the velocity profile, and calculate a local 
viscosity from the ratio of the two, it follows that this lo- 
cal viscosity is small near the moving inner cylinder, and 



increases with increasing distance from the moving wall. 
This is in fact easy to imagine, if the particle concentra- 
tion is slightly smaller near the wall, and increases with 
the radius. The expected concentration profile can easily 
be evaluated with a Krieger-Dougherty-like model for the 
dependence of the viscosity on the volume fraction; it is 
outside of the scope of this paper, but we plan to mea- 
sure the concentration profile with the MRI and compare 
them to what is expected from the velocity profiles. 

In conclusion, we have studied the transition between 
the frictional and lubricated flows of a dense paste. 
We find that the bifurcation of viscosity completely 
and unambiguously characterizes the transition between 
jammed and flowing states under imposed stress, and 
between localized flows and homogeneous flows without 
localization under imposed shear rate. The bifurcation 
gives both the critical stress and critical shear rate, the 
critical shear rate being inversely proportional to the fluid 
viscosity. We have shown that the critical Leighton num- 
ber, which follows from the critical shear rate and stress 
and from the solvent viscosity, completely fixes the tran- 
sition between the two regimes. 

In the 'frictional' regime, the shear flow is localized, 
with a shear zone increasing with the macroscopic shear 
rate unlike what happens for dry granular flows. More- 
over, the 'lubricated' regime is not simply viscous: the 
ratio of normal to shear stresses remains constant, the ve- 
locity profiles are roughly exponential and can be rescaled 
in a universal way in the two regimes. Finally we show 
that either the material is inhomogeneous, or that no sim- 
ple constitutive relation exists in the lubricated regime; it 
is in any case impossible to define a macroscopic viscos- 
ity that does not depend on the measurement system or 
the geometry, both in the frictional and in the lubricated 
regime. 
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